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The two-fold singularity has played a significant role in our understanding of uniqueness 
and stability in piecewise smooth dynamical systems. When a vector field is discontinuous 
at some hypersurface, it can become tangent to that surface from one side or the other, 
and tangency from both sides creates a two-fold singularity. The flow this creates bears a 
superficial resemblance to so-called folded singularities in (smooth) slow-fast systems, which 
arise at the intersection of attractive and repelling branches of slow invariant manifolds, 
important in the local study of canards and mixed mode oscillations. Here we show that 
these two singularities are intimately related. When the discontinuity in a piecewise smooth 
system is blown up or smoothed out at a two-fold singularity, the resulting system can be 
mapped onto a folded singularity. The result is not obvious, however, since it requires the 
presence of nonlinear or ‘hidden’ terms at the discontinuity, which turn out to be necessary 
for structural stability of the blow up (or smoothing), and necessary for mapping to the 
folded singularity. 


I. INTRODUCTION 


For a local singularity that is so easy to define, the ‘two-fold’ singularity of piecewise-smooth 
dynamical systems has proven surprisingly difficult to characterize. It has been the subject of 


interest both for its intricate phase portrait its ambiguous stability [10|, [ill], and recently 

because of its role in determinacy-breaking [2]. Now that it seems these are all well understood, 
our interest in this paper is in how the local dynamics relates to that of smooth flows. 

If a discontinuity occurs in the vector field of a flow on some hypersurface S, and the flow curves 
(or ‘folds’ parabolically) towards or away from S on both sides of the surface, then under generic 


conditions the resu 
Starting with 


t is a two-fold singularity, as depicted in figure [H 

111 ] , the elucidation of the two-fold singularity’s dynamics can be traced through 
references in (where the local form of the two-fold singularity is also extended to more than three 
dimensions). Though it can create limit cycles via a local bifurcation Q], the two-fold singularity 
itself is not an attractor. The main interest therefore lies in its transitory effects on a flow, 
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FIG. 1: Three kinds of two-fold. The main figure shows the phase portrait: in the unshaded regions the flow 
crosses (cr.) through a discontinuity at xi = 0, in the shaded regions the flow can only slide along the discontinuity 
on Xi = 0, the region being attracting {a.si.) or repelling {r.sL). In the examples shown, determinacy-breaking 
occurs at the singularity, meaning that the flow there becomes set-valued, the set has 2 dimensions in (i) and 3 
dimensions in (ii-iii). 


particularly in cases where it forms a bridge from an attracting region on the discontinuity surface 


into a repelling region, mimicking canard behaviour of smooth two-timescale systems [l| but with 
a more extreme character, because the two-fold breaks determinacy in both forward and backward 
time through the singularity. This is illustrated in figured] by a typical single trajectory entering 
the singularity, being deterministic until it does so, and afterwards exploding into a set-valued flow 
of infinite onward trajectories whose local form, however, is still somewhat constrained by the local 
geometry. 

Particularly because of some similarity to canard dynamics, attention has turned to how the two¬ 
fold can be understood as a limit or approximation of a smooth flow. The topological equivalence 
between “sliding” motion along a discontinuity surf ace, and “slow” motion on invariant manifolds 
of smooth two-timescale systems, has been shown [13l |. A qualitative connection has also been 
made between the sliding phase portraits at two-fold singularities and the slow manifold dynamics 
at so-called ‘folded’ singularities of two-timescale systems % ISJ. We shall prove here a more 
direct connection between the two singularities, by showing equivalence under explicit coordinate 
transformations. 

The two singularities that concern us can be defined as isolated points in n > 3 dimensions 
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satisfying the following systems of equations: 



The two-fold singularity: 

The folded singularity: 


system : 

/ 

Xi = /i(xi,X2,...;A) 

< 

for i = 1, ..., n 

< 

£ 2/1 = 51 ( 2 / 1 , 2 / 2 , •■•;£) 

Vj = 5j(2/i, 2/2, •••;£) 
for j = 2,..., n 

(1) 

parameter : 

A = sign (xi) 


0 < e < 1 


singularity : 

0 = xi = lim xi 


0 = 5i = 5i = 



along with a number of non-degeneracy conditions, which state that immediate higher derivatives 
have typical (i.e. non-vanishing) values; these will be elaborated on later. The functions fi and 
gi are assumed to be differentiable with respect to their arguments. The parameter A introduces 
a discontinuity via the sign function, while e introduces a separation into slow and fast timescales 
t and t/e (where i of course denotes the time derivative ^z). Three dimensions are sufficient to 
understand the basic dynamics of either singularity, so henceforth we restrict to n = 3. 

In an attempt to relate the two systems in ([T|) , we might first observe that the folded singularity 
system becomes discontinuous in the singular e ^ 0 limit. However, that limit is not equivalent 
to the discontinuous system of the two-fold singularity in ([T]) . To understand how the two systems 
in ([T|) are related we must instead start with the two-fold system, and somehow regularize the 
discontinuity to understand what forms of smooth system have it as a limit. 

Regularizing or smoothing a discontinuity raises immediate issues of uniqueness, namely that 
infinitely many qualitatively different smooth systems can have the same discontinuous limit. We 
shall find that actually this ambiguity can be encoded into the discontinuous system by means of 
nonlinear switching terms — nonlinear dependence on the parameter A in ([1]) — as introduced in 
P] (sometimes called ‘hidden’ terms because they vanish everywhere except at the discontinuity). 
We shall see that the smooth system obtained from a two-fold singularity is structurally unstable 
if it depends only linearly on A, and we shall show that a small perturbation, by terms that 
are nonlinear in A, restores structural stability and allows transformation into the general local 
singularity expected in a smooth system, namely the folded singularity. 

In section |II] we introduce the normal form of the two-fold singularity, and outline the basic steps 
for its study by blowing up the discontinuity in section IIIIl In section IIVI we blow up the normal 
piecewise smooth system, assuming only linear dependence on A as in the standard literature, giving 
a system which we find is structurally unstable. In section |V] we introduce a nonlinear perturbation 
of the normal form, then blow this up to find that it removes the structurally instability, and can be 
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mapped onto the folded singularity of a smooth two timescale system. Remarks showing that these 
results follow also if we smooth, rather than blow up, the discontinuity, are given in section IVTl 


II. THE TWO-FOLD SINGULARITY 


The normal form of the two-fold singularity is 


(xi,X2,i:3) = 


in terms of constants Oj = ±1 and bi G M. By results in 


(-X2,ai,5i) if xi > 0 , 

(+3^3, ^2, 02 ) if < 0 , 

flE, 


( 2 ) 


12l | , a system is locally approximated 


by ([2]) when it satishes the conditions in the lefthand column of ([T]), along with non-degeneracy 
conditions stating that the flow curvature is quadratic, given by lim ^ 0, and that the 

flows either side of the discontinuity are independent of each other and of the discontinuity surface, 
meaning the vectors along the xi direction, along lirn and along lirn xi\^^^_^, are linearly 


5^0 


independent. 

The local flow ‘folds’ towards or away from the switching surface xi = 0, along the line X 2 = 
xi = 0 on one side of the surface, and along the line X 3 = xi = 0 on the other. Hence the point 
where these lines cross is called the ‘twofold’. As a result, the surface xi = 0 is attractive in 
X2,X3 > 0 and repulsive in X2,®3 < 0 , while trajectories cross the surface transversely in X2®3 < 0 . 
In the attractive and repulsive regions the flow slides along the surface xi = 0, and the vector field 
it follows is found by interpolating across the discontinuity; we shall see how this is done below. 

The qualitative picture is then as shown in figure [TJ The precise form of the local dynamics 
depends on whether the flow curves towards or away from the discontinuity, determined by ai and 
02 , and also depends crucially on the quantity 6162 , which quantifies the jump in the angle of the 
flow across the discontinuity. An accounting of the many classes of dynamics that arise from these 
simple conditions is given in and references therein, we give only the pertinent details here. 

The three main ‘flavours’ of two-fold are: the visible two-fold for oi = 02 = — 1, the invisible 
two-fold for oi = 02 = 1 , and the mixed two-fold for 0302 = — 1 ; an example of each is shown 
in figure [H (i,ii,iii) respectively. In the cases depicted, there exist one or more trajectories passing 
from the attractive sliding region to the repelling sliding region, called canard orbits. This passage 
occurs in finite time (since the vector field ([2|) is non-vanishing everywhere locally). The flow is 
unique in forward time everywhere except in the repelling sliding region, where it is set-valued 
because trajectories may slide along xi = 0, but also be ejected into xi > 0 or xi < 0 at any point. 
This means that the flow may evolve deterministically until it arrives at the singularity by means 
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of a canard, at which point it becomes set-valued, so we say that determinacy breaking occurs at 
the singularity whenever canards are exhibited. This occurs in the invisible case when 61,62 < 0 
and 6162 > 1, in the visible case when 61 < 0 or 62 < 0 or 6162 < 1, and finally in the mixed case 
when 61 < 0 < 62 and 6162 < —1 or when 61 -|- 62 < 0 and 61 — 62 < —2. (The particular cases 
shown in figure [H are: (i) ai = 02 = —1 with 61 < 0 or 62 < 0 or 6162 < 1; (ii) oi = 02 = 1 with 
61,62 < 0 and 6162 > 1; (hi) 0102 = —1 with 61 < 0 < 62 and 6162 < —1 or with 61 -|- 62 < 0 and 
61 — 62 < —2.) 

We shall henceforth be concerned only with the dynamics that has previously been missing, that 
which occurs inside the discontinuity itself and gives rise to determinacy-breaking at the singularity. 
In any event, the first step in studying the local dynamics must be to extend the system ([2]) in 
such a way that is can be solved across the discontinuity xi = 0, as follows. 


III. BLOW UP 


We shall outline the procedure for studying ([2|) in a general form first. Let x = (xi,X2,X3) and 
f = ifi, f 2 , fa), and begin with a piecewise-smooth system 


X = 


f(x; -1-1) if xi > 0 , 
f(x; —1) if xi < 0 , 

which can be extended to xi = 0 by defining f as 

f(x; A) = ^-^f(x; +1) + i^f(x; -1) + {1- A2)g(x; A) 


(3) 


(4) 


where 


A € 


sign (xi) if xi 7^ 0 , 


(5) 


—1, -|-1] if xi = 0 . 

The function g is smooth, and provides f with a nonlinear dependence on A. Without g the 
unction f(x;A) is the convex combination as used in Filippov’s standard approach to switching 
8|, and here we shall introduce g 7^ 0 only as a small perturbation when necessary. 

The dynamics outside the discontinuity (xi 7^ 0) is now given by taking either of A = ±1 in Q, 
which reproduces ([3|). On the discontinuity (xi = 0) we consider whether there exists A G [—1, -|-1] 
such that xi = 0, for which the flow of (jH) lies in the tangent plane of the discontinuity, in which 
case sliding dynamics along the discontinuity surface is possible. The sliding dynamics can be 
reached in finite forward or backward time, meaning the flow can lose or gain a dimension as it 
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arrives at or departs from xi = 0, so the resulting flow might not be unique. Therefore we shall 
follow Q] and regain uniqueness by blowing up the discontinuity manifold xi = 0 , to study the 
dynamics of A that transports the flow through the jump between —1 and +1. This is done as 
follows. 

Because A is a function of xi only, the dynamics of A is induced by the xi component of the 
flow such that A' = /i(0, X 2 , X 3 ; A), where the prime denotes differentiation with respect to some 
dummy timescale, instantaneous compared to the main timescale. One way to describe this is to 
say A' = for infinitesimal e > 0. (In fact, this blow up is equivalent to defining A as a function 
A = 4>{v) of a fast variable v = x\/e such that (^[xxje) —>■ sign(xi) as e ^ 0; we discuss this further 
in section EH). We obtain for the blow up system on xi = 0, 

(A',i: 2 ,X 3 ) = (/i(0,X2,X3; A), /2(0, X 2 , X 3 ; A), /3(0, X 2 , X 3 ; A)) . ( 6 ) 

The steady states of the ‘fast’ prime system satisfy a differential-algebraic system 

(0,X2,X3) = (/i(0,X2,X3; A), /2(0,X2,X3; A), h{0,X2,X3;X)) , (7) 

which defines dynamics on an invariant manifold 

= {(A,X2,X3) G X : /i(0,X2,X3; A) = 0} . 

The system ([7]) thus defines sliding dynamics in the discontinuity surface xi = 0, fixing the value 
of A G [—1, -t-1] inside the manifold (when it exists) for which xi = 0 has invariant dynamics, 
and specifying the variation (x 2 ,X 3 ) on that manifold. In regions of xi = 0 where does not 
exist (where /i / 0 for xi = 0 and A G [—1,+!]), the system ([ 6 |) conveys the flow across the 
discontinuity from one region xi ^ 0 to the other. 

The systems dH and ([ 6 ]) are sufficient to specify the local dynamics, and we shall now apply 
these to the normal form. The final step in our analysis will be to transform the system ([ 6 ]) into 
the normal form of a folded singularity, but we shall find that this is only possible if g is nonzero, 
and hence there is a nonlinear dependence on A. 

IV. THE UNPERTURBED SYSTEM 

In this section we show the following. 

Proposition 1. The blow-up ^ of the normal form two-fold singularity ^ with g = 0 is struc¬ 
turally unstable. 
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To prove this we will perform the blow up described in section [H and show that the sliding 
manifold exhibits a degeneracy. 

To resolve the dynamics across the discontinuity, we use (jl]) to express the system ([2]) as a 
convex combination 

{ xi , X 2 , X 3 ) = ^-^(-X2,ai,6i) + ^-^(0:3,62,02) (8) 

:= (/i(xi,X 2 ,X 3 ; A), / 2 (a:i, X 2 , X 3 ; A), / 3 ( 2 :i, X 2 , X 3 ; A)) . 

Using the blow-up ([ 6 ]) on xi = 0 we obtain the two timescale system 

(y,X 2 ,X 3 ) = (/i,/2,/3) = ^^^(- 3 ^ 2 , 01 , 61 ) -h ( 3 : 3 , 62 , 02 ) , (9) 

with A € [—1,+!]. The fast system on A^ has equilibria where A^ = 0, yielding a differential- 
algebraic system 

(0,X2,X3) = (/l,/2,/3) = ^-^(-X 2 ,ai, 6 i) -h ^-y^(x3,62,a2) . (10) 

This describes states that evolve inside the switching surface xi = 0 on the main timescale, and 
these are precisely Filippov’s sliding modes. We can solve in the sliding mode to get A = (xs — 
3 : 2 ) 7 ( 3:3 + X 2 ), and sliding modes exist only when this lies in the allowed range A G [—1,-|-1], hence 
for X 2 X 3 > 0 . 

In the absence of sliding modes, where X 2 X 3 < 0, equation m describes the instantaneous 
transition from one boundary of A G [—1, -|-1] to another, whereby the flow crosses through the 
switching surface as A flips sign. 

Concerning the sliding regions X 2 X 3 > 0 on xi = 0, the subsystem m inhabits invariant 
manifolds of the blow-up system ([9]) given by 

= \ (A, X 2 , X 3 ) G [-1, +1] X : A = ——xi = 0 < X 2 X 3 1 . (11) 

I X 3 + X 2 j 

We call the sliding manifold. It consists of two normally hyperbolic branches, one attractive 

in X 2 ,X 3 > 0 since dX'/dX = —{X 3 -|- X 2)/2 < 0, and one repulsive in X 2 ,X 3 < 0 since dX'/dX = 

— (x 3 -|- X 2)/2 > 0. The two branches are connected at X 2 = X 3 = 0, along a line 

C = {(A,X 2 ,X 3 ) G [-1,+1] X : X 2 = X 3 = 0, } . (12) 

Along C we have dX'/dX = 0, so the manifold is no longer normally hyperbolic where it meets 
C, and therefore is not invariant there. Figure [2] shows an example of the discontinuous system (i), 
and its blow-up (ii) showing Ai^ and £, rotated in (iii) to show £ more clearly. 
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regularize 


rotate 




FIG. 2: Blowing up the unperturbed system for the example of an invisible two-fold, (i) The flow directions 
outside a;i = 0 create an attracting sliding region in X 2 ,X 3 , > 0 and repelling sliding region in X 2 ,X 3 < 0. (ii) The 
blow up of xi = 0 into the interval Ai G [—1,-|-1], where the sliding regions become an invariant manifold 
(shaded), hyperbolic except along the vertical line C, which aligns with the fast (double arrowed) Ai dynamics, (iii) 
The dynamics in the manifold is best viewed along the U3 axis of rotated coordinates U2 = X2+X3, U3 = X2 — X3. 


The line L represents the two-fold singularity (xi = X2 = X3 = 0 ), stretched out over the interval 
A G [— 1 ,-|- 1 ] in the blow up. Comparing the alignment of L with the dynamics outside , 
however, reveals a degeneracy in the system ([ 9 |). The non-hyperbolic set £ is a line with tangent 
vector ec = (1,0,0) in the space of (A,X2,X3), which means it lies everywhere parallel to the fast 
one-dimensional A^ system in ([ 9 ]). As a result, all along £ we have not only dfi/dX = 0 , but also 
d'^fi/dX^ = 0 and moreover d^fi/dX^ = 0 for all r > 1 . This constitutes an infinite codimension 
degeneracy, so the system is structurally unstable. This is the result in the proposition. 


Remarks on the degeneracy 

The degeneracy means that the flow is not transverse to the set £ where the two sheets of 
intersect. As a result, if we perturb the system, transversal intersection between the attracting and 
repelling branches of the invariant manifold cannot be guaranteed, see e.g. 141 . Small changes 
in the expression ([8]) may result in a system where the two sheets of intersect at points, along 
intervals, or where no intersections occur at all, dependent on the kind of perturbation. These 
intersections are particularly important because they are they points where canard trajectories can 
flow from the attracting to repelling branches of Ad'^ (or faux-canard trajectories can flow in the 
opposite direction). Hence in the system above the existence of canards cannot be guaranteed. 

In the literature on smooth two timescale systems (as we have here in (A,X2,®3) space), the 
connection of attracting and repelling branches of a slow invariant manifold has been well studied, 
leading to a generic canonical form and requisite non-degeneracy conditions as described in Q- 
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In our notation those conditions are the existence of a point on £ where 


/i = 0 


dfi _ n 


dfi 

dX2,3 


/O 


( 13 ) 


the fourth of which is violated everywhere in the system described above. 

Adding constant terms or powers of X2 and X3 to ([8]) would only move the set £ in the (x2,X3) 
plane, not remove its degeneracy, and terms depending on xi would vanish at the discontinuity. 
The only recourse to break the degeneracy, particularly to give / 0 , is therefore to add terms 
nonlinear in A. Anything we add to the function fi in ([8]) must still give ([2|), so it must vanish 
outside the switching surface xi = 0, i.e. be a perturbation of the form dH), as we introduce in the 
next section. 


V. THE PERTURBED SYSTEM 

In this section we show the following. 

Proposition 2. The blow-up ^ of the normal form two-fold singularity ^ with g = (a, 0 , 0 ) can 
be transformed into 

exi = X2 + x? + 0 (exi,ex3,xix3) 

X2 = bX3 + CXi + 0 (x|, X1X3) 

X3 = a + 0 (x3,xi) 

provided a / 0 for small e > 0 , where a, 6, c, are real constants, and provided the conditions 
^(61 — 62) < 1 = ai = —a-2 or ^{bi — 62) > —1 = oi = —02 do not hold. 

It turns out that the case excluded by the conditions ±^(61 — 62) < 1 = ±ai = =F«2 is that in 
which there are no canards or faux-canards, i.e. no orbits of the sliding flow passing through the 
singularity. 

The perturbed system we consider, applying ([ 1 ]) to the normal form ([ 2 |) with g = (a, 0 , 0 ), is 
(xi,X2,®3) = {-X 2 , ai,bi)+ ^ (x3,62, 02) + (1 - A^)(a, 0, 0) ( 14 ) 

:= (Fi(xi,X 2,X3; A),F2(a;i,X2,X3;A),F3(xi,X2,X3;A)) , 

where a is a constant. We shall show that perturbing xi with a small term a in this way is neces¬ 
sary and sufficient for structural stability. Perturbing X2 or X3 is neither necessary nor sufficient. 
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therefore we leave them unaltered with ^2,3 = /2,3 (referring to {fi, f2, fs) in ([S])). The systems 
(IH), ([8]) and (fTT|l are identical for all a: 1 7^ 0, when A = ±1. 

On = 0 we consider the blow up of (I 14 jl as before, now giving 


(A', X2, X3) = i-X2, ai, bi) + (x 3 , 62,02) + ot{l - A^) ( 1 , 0 , 0 ) 

This is an a-perturbation of Q. The system has equilibria on the set 


— < (A,X 2 ,X 3 ) G [— 1 ,+ 1 ] X 


1 + A 1 -A , 

-X2 H-—X3 + a{l - A ) = 0 


( 15 ) 


( 16 ) 


2 " 2 

The surface described by forms an invariant manifold of (| 15 l) at all points where it is normally 


o \ / 

hyperbolic, i.e. where ^ 7^ 0, which now holds except on the curve 


C = < {\,X2,X3) : A = 2 


2a + X3 - X2 X3 + X2 


G [-1,+1] 


( 17 ) 


X3 + X2 4 a 

The curve C has tangent vector e/; = (1 , + 2 a(A — 1 ) , — 2 a(A + l)), which for all |A| < 1 is 
transverse to the flow if a 7^ 0 . In particular L is now transverse to the dot and prime timescale 
snbsystems in (| 15 p . so the degeneracy of the unperturbed system has been broken. We show below 
that this perturbed a 7^ 0 system is structurally stable by mapping the two timescale blow-up 
system m onto the generic folded singularity of smooth systems. 

While the non-hyperbolic curve C is now in a general position with respect to the flow, there 
may exist a new singularity along £, where the flow’s projection along the A-direction onto the 
nullcline /i = 0 is indeterminate. This is the so-called folded singularity, found as follows. The set 
/i = 0 is a fixed point of the A' subsystem of (I 15 p the flow. A trajectory which remains on /i = 0 for 
an interval of time satisfies 0 = /i = (A', X2,X3) ■ {d/dX, d/dx2,8/8x3) fi = + (/2, fs) • 

which rearranges to ^y-A' = —(/2, fs) ■ 3:3) ' indeterminate at points where 

dfi 


Denoting the value of fi at the singularity as fis, and solving (fT8l) via 

1 + As 1 — As 


( 18 ) 


0 — {f2s, hs) ■ — 


— _|_ ai2lh.X^ fei+Qa _j_ b-i -a-2 


lH~As 1—Ag 
2 ’ 2 


( 19 ) 


we find that the folded singularity lies at {X,X2,X3) = (As,X2s,a^3s) where 

ai+a2 /i I 4aia2 

^ _ 61-62 V _ _ I „Y\ 1^2 „ _ , ii 2 

As — , ^2s — Ta(As 1) , X3s — a(As -I- 1) 


1 I Ql —Q2 
^ 61-62 


(20) 


Noting that oi and 02 jnst takes values ± 1 , we have: 
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in the case ai = 02 = 1, = 5^-:^ =t there exists a unique A* G [—1, +1] 

for any bi and 62 (the positive root for bi >62, the negative root for bi < 62); 


in the case oi = 02 = — 1 , A* = ± ^1 + , so there exists a unique As G [— 1 , + 1 ] 

for any 61 and 62 (the positive root for bi <621 the negative root for bi > 62); 


• in the case oi = —02 = 1 , As = there exist two points with As G [— 1 , + 1 ] for 

61 — 62 > 2, and no points otherwise. 

• in the case ai = —02 = — 1 , As = ± there exist two points with As G [— 1 ,+ 1 ] for 

61 — 62 < —2, and no points otherwise. 

In the cases where As is unique we proceed directly to the steps that follow below. In the cases where 
As can take two values we can proceed with the following analysis about each value, and will obtain 
different constants in the final local expression, i.e. a different folded singularity corresponding to 
each As. In the cases when As does not exist, no equivalence can be formed; these are the cases 
when the the two-fold’s sliding portrait is of focal type (see j^) where orbits winds around the 
two-fold but sliding orbits never enter or leave it. So excluding those cases ai = —02 = 1 with 
61 — 62 < 2 and ai = —02 = —1 with 5 i — 62 > — 2 , we can proceed as follows. 

Taking a valid value for As, a translation puts the singularity at the origin of the coordinates 


then fi becomes 


2/1 — A — As , y 2 — X 2 — X 2 s , Us — X 3 — XSs , 


As -k 1 As — 1 / 2/3 + 2/2 , , 

Ji =-—2/2-—2/3 - --k ayi yi 


2 2 

with derivatives evaluated at the singularity 

^ _ _m±m _ 2(ym 
dyi ~ 2 ’ 

S# = 


( 21 ) 


( 22 ) 


2//1 _ _ f Ag-l-l As —1 1 

~ V 2 ’ 2 / ’ 


diy2,y3) V 2 ’ 2 

^ 


( 23 ) 


d{y2,y3)dyi 

Solving the conditions fi = = 0 we find that C can be written in functional form as 

( 2 / 1 ,2/2, 2 / 3 ) = ^( 2 / 1 ) := (2/i,-a2/i(2-2As-2/i), -ayi(2-k 2As-k 2 / 1 )) • 


It is more useful to parameterize £ in terms of ys, which we obtain as 

^ yiiiys) \ ( -1 - As -k v^(l -k As)2 - ya/a ^ 


C : 


2 / 2 l( 2 / 3 ) 

V 2 / 3 l( 2 / 3 ) y 


- 2/3 - 4 a(-l - As -k y/(l -k As)2 - ys/a) 
2/3 


( 24 ) 
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in terms of three functions yiL,y2L,y3L, whose derivatives are 

> f ^ -l/ 2 a , , X _ 1 - As- 2 /iz,(y 3 ) , 

- TTArr^Hfe) ■ + + 

We can then rectify C to lie along some Z3 axis by defining new coordinates 
zi = yi- yiiiys) , ^2 = 2/2 - y2L(y3) , Z3 = ys . 


( 25 ) 


( 26 ) 


The former vector field components become 

^ _ i+K+yiL ^^ _ ^^2 _ , 

/2 = a+k + a^As + 2i^(zi+yiL(z3)) = /2s + 0(z3,^i) , (27) 

/s = 2^^ As + 2 ^^ (^1 + yii(^3)) = /ss + 0 (Z3, Zi) . 

When we now evaluate the derivatives z[,Z2,Z2, it is necessary to relate the dot and prime time 
derivatives. We do this using an infinitesimal e > 0 , letting a dot denote the derivative with 

respect to t, while a prime denotes the derivative with respect to a ‘fast’ time t/s which becomes 

instantaneous as e —)■ 0 . Then 


4 


fi - £Z3y[Liz3) 

As +1 / e/3s \ 

“r'"a(l + As) 2 y' 


azf + 0{sZ3,eZi,Z2Z3,Z2Zi) , 


using yiL{z3) = 0 {Z3) 
which yields 


— + 0 (z|) . We then make a small corrective shift Z2 = Z2 — 


sfss 

a(l+Aa)^ 


^2 = /2 - hy'2L{z3) = -Z3 + CZi + 0 (zf, 2 ; 32 ;i) 

a 


where 

df2s 1 - As df3s 1 / /2s + hs \ 

dX 1 + As 5 A ’ “ 2 V 1 + As ^ 7 ■ 

The last thing to do is just scaling. Collecting everything together so far we have 


z[ = diZ 2 - azf + 0 {ez,ez 3 ,ziZ 3 ) 

Z 2 = ^Z3 + CZI + 0 (zf, 2 ;i 2 ; 3 ) ( 28 ) 

7 = fss + 0(2:3, zi) 

where di = —^(1 + As). Defining new variables xi = X2 = —s\gii{a)diZ2-i X3 = —sign(Q;)2;3, 

i = —sign(a)t, gives 


x[ = X2 + xl + 0 {exi,ex3,xix3) 
X2 = bx3 + cxi + 0 [xl,XiX3) 

X 3 = a + 0(x3,xi) 


(29) 
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where 

a = fss , b = (f2s + fss - 2cv^) , 

af \ ( 30 ) 

+ + • 

Since x' = £x this is the result in the proposition, and it is clearly valid only for a 7^ 0 . 

Figure [ 5 ] shows an example of the perturbed system and its blow up for each flavour of two-fold 
in (i) (corresponding to those in figured]), followed by their blow up (ii), and a rotation (hi) to 
show the phase portrait around the set £ more clearly (similar to figure [ 2 |). In the most extreme 
case, the folded node, the original phase portrait contains infinitely many intersecting trajectories 
traversing the singularity, while the perturbed system splits these into distinguishable orbits, a 
finite number of which asymptote to the attracting and repelling branches of the critical manifold. 


Remarks on the singularity 


A glance at the papers Q, 14 , Q reveals what a charismatic singularity lies hidden in the 
dynamics of the two-fold, revealed by blow-up, and waiting to be released when the discontin¬ 
uous system is perturbed by simulations that smooth, regularize, or otherwise approximate the 
discontinuity. 

In particular one may ask what happens in the cases of determinacy-breaking illustrated in 
figure dl since a smooth flow should be deterministic. Because the two-fold maps to a folded 
singularity (actually one or two folded singularities as we saw above), if we consider what happens 
with e > 0 we find determinacy is restored, but replaced by sensitivity to initial conditions. 

The manifold becomes the slow critical manifold of a slow-fast system, and for e by standard 


0 , 


there exist invariant manifolds A 4 f in the 


results of geometric singular perturbation theory 
e-neighbour hood of but away from C, on which the dynamics is topologically equivalent to 

the sliding dynamics (termed ‘reduced’ or slow dynamics usually in this context) found above. 
Trajectories that pass close to the singularity, or more precisely, close to the folded singularity 
on the set £, may persist in following the manifold from its stable to unstable branches, 
while other nearby trajectories will veer wildly away, their fate sensitive to initial conditions and 
proximity to primary canard orbits (those which persist along both branches of throughout 
the local region). 

Like the different kinds of two-fold, there are different classes of folded singularity, and their 
classification depends on the slow dynamics inside A 4 ^. From the expressions (| 29 p - () 30 p we see that 
the class therefore depends not only on the constants ai, 02, j ^2, of the original piecewise smooth 
system, but also on the ‘hidden’ parameter a. 
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FIG. 3: Blowing up the perturbed {a ^ 0) system, for examples of each flavour of two-fold. Labelling as in 
figure [2] Note in the blow up (ii) that C is now a curve. Rotating around the A axis in (iii) we can see the 
attracting branch (upper right segment) and repelling branch (lower left segment) of the invariant manifold 
(shaded), connected by £. The folded singularity (f.sing.) appears along £, two in the case of mixed visibility, 
recognised as having a phase portrait that resembles a saddle or node if we reverse time in the repelling branch 
of . In (iv) we sketch the corresponding phase portraits in the slow-fast system ()29p . 


The classification scheme is fairly simple, and can be used to verify the dynamics on seen 
in figure [ 3 l The projection of the system (j 29 p onto , found by differentiating the condition 
f) = X2 -\- x\ with respect to time to give 0 = 6x3 -|- cxi -|- 2 xixi -)- 0 (x|, X1X3), is 


Xl 

^3 


1 

—2xi 



-h 0 (x|,XiX3) 


A classification then follows by neglecting the singular prefactor l/ 2 xi and considering whether 
the phase portrait is that of a focus, a node, or a saddle. This is determined by the 2 x 2 matrix 
Jacobian, which has trace c, determinant 2 ab, and eigenvalues 5(0 ± ■\/c^ — 8 db). This will not be 
the true system’s phase portrait because the time-scaling from the l/2xi factor is positive in the 
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attractive branch of A 4 ^, negative (time-reversing) in the repulsive branch, and divergent at the 
singularity (turning infinite time convergence to the singularity into finite time passage through 
the singularity). The effect of this is to ‘fold’ together attracting and repelling pairs of each 
equilibrium type, so each equilibrium becomes a ‘folded-equilibrium’, forming a continuous bridge 
between branches of . As a result the flow on is a folded-saddle if a6 < 0, a folded-node 
if 0 < 8o6 < (?, and a folded-focus if (? < 8ab. Canard cases occur for c > 0 and faux canard for 
c < 0 . In figure [ 3 ] we show the result of blowing up the discontinuous system for an example of each 
type of two-fold that exhibits determinacy-breaking (those from figured]). In the visible two-fold 
the singularity becomes a folded-saddle, in the invisible case it becomes a folded-node, while the 
mixed case becomes a pair consisting of one folded-saddle and one folded-node. 

One may ask why certain cases, in particular those without determinacy-breaking, were excluded 
by the proposition above. These are systems where the sliding phase portraits contain no orbits 
passing between the attracting and repelling sliding regions, and therefore one expects that in 
their blow up there should exists no canards connecting attracting and repelling branches of 
and hence no folded singularities. If there are no such connections then no intersections between 
the two branches of persist when we let e be nonzero. Hence the omission of these cases is 
consistent, and a posteori it is obviously necessary. 

In figure 0 ] we illustrate the way splitting of manifolds occurs when we introduce the instanta¬ 
neous timescale (whose derivative is denoted with a prime) with a fast timescale t/e for e > 0 . The 
case shown is an invisible two-fold. 



FIG. 4: Sketch of: (i) a two-fold singularity; (ii) its blow up with o yf 0; and (iii) the splitting of the branches 
of when we replace x[ with exi and consider small e > 0. (iii) approaches (ii) as e — >■ 0. 

An in-depth description of the dynamics that ensues in the different cases of two-folds, and 
the smoothings subject to perturbations, would be lengthy, and deserves future study elsewhere. 
As a demonstration we conclude with three examples of complex oscillatory attractors formed by 
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two-fold singularities, simulated by smoothing. We take a = 1/5 and: 


(i) f+ = (-X2, fiCl -k ^X2 - 1, ^X2 - \x 2 X 2 , - |), f~ = (X3, \x2Xz - f, |x3 - 1 - Xi)] 

(ii) f+ = (-X2,1 -k Xi, -|), f“ = (X3, 1 - f Xi); 

(hi) f+ = (-X 2 -k ]^xi,xi - |,xi - 2 ), f” = (x 3 -k ^xi,xi -k 1 - xi); 


simulated in figure [ 5 j The simulation uses Mathematica’s NDSolve, for which the discontinuous 



FIG. 5: Three examples of attractor organised around a two-fold singularity. Showing: (i) a simulated trajectory, 
(ii) a sketch of the discontinuous flow inside and outside xi = 0 that gives rise to it, and (iii) the blow up on 

xi = 0 . 

function Ai = sign xi is replaced with a sigmoid function tanh(xi/e) with e = 10 “®. Because the 
system is structurally stable, the exact choice of smoothing function is not important so long as 
it is simple, for example continuous and monotonic, and further simulations not shown here verify 
that such dynamics persists with different smoothing functions, for example varying e, or replacing 
xi with a different sigmoid function such as xij+ xf. The hrst row shows the simulation of a 
single trajectory for a time interval t = 1000 in (xi,X2,X3) space, as the flow attempts to switch 
between the critical slow manifold flow in (ii) given by the blow up on xi = 0, and the xi > 0 and 
xi < 0 flows in (i). In the second row we plot the familiar picture of sliding dynamics on xi = 0 
for the unperturbed system (neglecting nonlinear A), and sketch the curvature of the xi 7^ 0 flows. 
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The sliding vector field has canard trajectories passing from the righthand attracting branch to the 
lefthand repelling branch via the folded singularity; at a visible two-fold only one canard exists, 
at an invisible two-fold every sliding trajectory is a canard, and at a mixed two-fold a region of 
trajectories are canards. The third row shows the critical manifold in the (Ai, X2, X3) blow up. 


VI. CLOSING REMARKS 


For such a simple system, even taking its piecewise linear local normal form, the two-fold 
exhibits intricate and varied dynamics. Just how intricate becomes even more clear as we attempt 
to regularize the discontinuity, and study how the two-fold related to slow-fast dynamics of smooth 
systems. As well as insight into the dynamics that is seen upon simulated such a system, this 
adds a new facet to the question of the structural stability of the two-fold, which has remained a 


stimulating question since 
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The degeneracy in section ITVl gives some insight into why simulations of systems containing two¬ 
fold singularities with determinacy-breaking singularities are highly unpredictable. In attempting 
to solve near a discontinuity, simulations may introduce aspects of smoothing, hysteresis, delay 
and noise, not represented in the ideal discontinuous model, with unpredictable effects on such a 
structurally unstable high codimension degeneracy. We therefore have to resolve this degeneracy 
before expecting to obtain comprehensible simulations, and without leaving the class of piecewise 
smooth systems, the nonlinear switching term introduced in the perturbed system in section |V] is 
sufficient. 

The system (flTl) now succeeds ([8]) as our prototype for the local dynamics of the two-fold singu¬ 
larity (with A given by ([ 5 ]) in both cases). The question of whether this constitutes a ‘normal form’ 
has issues both in the discontinuous and slow-fast settings, but it is clear that (H is structurally 
stable, and represents all classes of behaviour that occur both in the discontinuous system, and in 
its blow up to a slow-fast system. 

Finally we should emphasize that the results of this paper are equivalent if we smooth, rather 
than blow up, the discontinuity. In particular the degeneracy that makes the perturbation g 7 ^ 0 
necessary is still present. To relate the procedure of smoothing (or Sotomayor-Teixeira’s regular¬ 
ization [1^) to the blow up procedure, we consider A as a function rather than a dynamic variable. 
To smooth ([8]) we let A = 4 >{xi/e) for small e > 0, then define a new variable v = exi for which 
v' = £v = /i(eu, X2, X3; 0 (u)), which replaces A' in ([ 9 ]) to give us our two timescale system. Note 
that and v' are related by scalar, A' = eA = £(f)'{v)v, amounting only to a time rescaling. Hence 
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the (/>-smoothed system is qualitatively similar to the A-blow up system, with the same invariant 
manifold A4^, and containing the same line L which lies aong the fast v' flow when e = 0. 
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